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Abstract 

We review the existing results on the scaling dimensions of oper- 
ators with more than two derivatives in the non-linear sigma models. 
We argue that the speculations on the relevance of these operators, 
and correspondingly on the breakdown of the (d— 2) expansion for the 
classical Heisenberg model, or for the one-parameter scaling theory of 
localization, are based on a dubious mathematical analysis. 



The low temperature expansion of the A- vector model is given by a model 
of interacting Goldstone bosons, known as the non-linear sigma model |], §]■ 
The existence of a zero-temperature fixed point in dimension two has led to an 
expansion of the critical properties in powers of (d— 2) which supplements the 
familiar expansion in powers of (4 — d) based on the linear Landau- Ginzburg- 
Wilson model 0]; linear and non-linear refer here to the representation of the 
0(A) symmetry by the order parameter. 

These non-linear sigma models were generalized to situations in which 
the order parameter is a mapping of a two-dimensional lattice onto a target 
manifold (a sphere in the case of the 0(A) model). These generalizations are 
the basis of the conformal field theory approach to string theory, and they 
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have been also applied to various problems in condensed matter physics, in 
particular to Anderson localization 0] , for which the manifold is the homoge- 
neous space 0(N, N)/0(N) x O(N) [§, §, 0], in the replica limit in which N 
goes to zero. This had led to an expansion of the metal-insulator transition 
in dimension three in powers of e = d — 2, with e = 1 at the end. The numer- 
ical success of this expansion has been unfortunately in no way comparable 
to that of the Wilson- Fisher expansion in powers of (4 — d). However for 
the localization problem, there is no upper critical dimension from which an 
expansion similar to the (4 — d) would allow us to ignore the poor numerical 
accuracy of the (d — 2) expansion. There is no other analytic approach to 
the problem. 

In view of these difficulties a growing number of articles, starting with the 
work of Kravtsov, Lerner and Yudson ||, extended later to the more familiar 
O(N) model by Wegner M, pointed out to the possible appearance of new 
relevant operators, which would of course drastically change the physics and 
the results of the model. 

In this note we would like to put the calculations which led to the alleged 
breakdown of the (d — 2) expansion, in the right perspective. Not that we 
have any significant additional calculation to report, but we would like the 
reader to judge the seriousness of the problem by presenting a summary of 
the situation as it is known at present. 

We shall center our discussion on the O(N) model whose physics is under 
much better control, and for which the 4 — d and 1/N expansion, plus a 
wealth of numerical experiments leave little doubts on its properties. Let us 
start by reminding the reader that the model is defined by a lattice of unit 
iV- vectors with nearest neighbour interactions of the form S(x) ■ S(x + ae a ) , 
in which e a is a unit vector in the direction a. In the long distance limit we 
can replace S(x) by a continuum field and the interaction by S- V 2 S, with the 
constraint S 2 = 1; (the term involving one derivative of S cancels with the 
x — ae a neighbour). The model is then defined by a straight low temperature 
expansion over some ordered state. However this expansion is plagued with 
short-distance divergences, as in any continuum field theory, and these diver- 
gences have to be regulated. This is slightly tricky in this problem; a straight 
ultra-violet cut-off would not work since it breaks the O(N) invariance of the 
model. There are two schemes known for regulating these divergences with- 
out breaking the rotational invariance, the dimensional regularization or a 
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lattice, which means simply returning to the low temperature expansion of 
the original classical Heisenberg model. Of course these two regularizations 
lead at the end to identical results in the critical long-distance limit In 
other words one should be clear at this stage: if one does find a problem 
with this theory it is not simply a problem of some abstract non-linear field 
theory, but of the low temperature expansion of the Heisenberg model itself. 

The difficulties which have bothered a number of workers concern the 
scaling behaviour of operators with more derivatives than two. Indeed it is 
clear from the very definition of the model that one could have operators 
such as S ■ (V 2 ) 2 >S, or more complicated ones which would involve derivatives 
and four spins for instance, as the ones which would be generated by the 
interaction of four spins around a plaquette. In any way we know from 
Wilson theory [|1(J that we should allow in the Hamiltonian for any possible 



O(N) invariant operator and find the fixed point in a space of an infinite 
number of coupling constants. Note however that we know a priori that 
there is one and only one relevant operator (in the absence of a symmetry- 
breaking operator such as a magnetic field); indeed it is sufficient, in order to 
reach the critical point to vary a single parameter, namely the temperature. 
If we had more relevant operators, we would have more "knobs" to adjust 
before we could reach criticality. Numerical experiments, expansions such as 
the 1/N, have never found anything of the sort: it is sufficient to be at T c to 
obtain an infinite correlation length. 

In two dimensions, simple dimensional counting of operators involving 2s 
derivatives of the order parameter ( an even number is needed to make a 
scalar), gives their scaling dimension 

y s = 2-2s (1) 

which leads to conclude that operators with more than two derivatives are 
irrelevant. There is no operator for s = since S 2 = 1; therefore we are left 
with a single relevant operator, namely S ■ V 2 S, as expected. 

In dimension d two things occur : (i) the canonical dimension is changed 
to (d — 2s) (ii) scaling anomalies appear, calculable from the renormalization 
group theory by expanding in powers of (d — 2). Wegner's result 0, based 
on a one-loop calculation, gave 

ys = d-2s+ € ^f^ + 0(e 2 ) (2) 
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and a recent two-loop calculation by Castilla and Chakravarty (TTJ gave for 
large s 

Vs = d-2s + ^fil + + 0(s 2 )] + 0( e 3 ) (3) 

(their calculation is done for finite s as well but the result is only quoted for 
large s). Finally let us quote the last available information which is the 1/N 
expansion of Vasil'ev and Stepanenko |12j which gives 



4 (d - l)T(d - 2)[s(s - l)d(d - 3) - 2s{d - 2)] 1 

Vs + n (4 - d)v{2 - f)r(f - i)r(f - 2)r(f + 1) + w 

i.e by expanding also in powers of e = d — 2 

We remark here that we could determine the two loop result of y s (||) 
explicitly by noting that y\ = 1/V, where v is a critical exponent for the 
correlation length, and that t/2 = —w, where u is an exponent for the cor- 
rection to the scaling. Indeed, (|J) for s = 2 agrees with the expression of 
1/N expansion of the exponent of the correction to the scaling by Ma [13 



From the results of 1/N expansion, and from (|j) the two- loop term may be 
expressed as (|s 3 + as 2 + (N + b)s)t 2 c . Since we know that it is proportional 
to (N — 2), for both s = 1 and 2, this fixes the two unknown constants a and 
b to be a = —2 and b = — |. 
Thus we get 

Vs = d - 2s + s(s - 1% + |s(a - - 2) + s{N - 2)]t 2 c + 0(t 3 ) (6) 

where t c is obtained from the zero of the (3 function [[J . Putting the expression 
for t c in the e expansion into (0), we have 

Vs = 2 + e-2s + s(s-V 



N -2 

2 , . w „ „ S1 e 2 



+ [-8(3 - 1)(5 - 2) - 8 ( 8 - 1) + 8 (iV - 2)] — _ + 0(6 3 ) (7) 
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We may recover this result by an alternative method: we consider again 
operators with 2s derivatives, but for sigma models on the target manifolds 
0(N)/0(p) x 0(N — p). From the isomorphism between 0(5)/ 0(4) and 
Sp(2)/Sp(l) x Sp(l) 0, we conclude that 

Vs = d _ 2s + s(s - l)t c + - 1) (s - 2) + (2p(N-p) - N)s]t 2 c + 0(t c 3 ) (8) 



using now 



(2p(N - p) - N)e 
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N-2 (N-2) 3 

we find 



+ 0(e 3 ) (9) 



y s = 2 + e-2s + s(s-l)j^-^ + [^s(s-l)(s-2) 

- MN ~*l~ N '(* ~ 1) + (MN-P) - N)s]y^yi + 0(e*) 

(10) 

which reduces to (0) for p = 1. 

The results of Wegner and followers led to question the possible relevance 
of the operators with more than two derivatives. Indeed, if in two dimensions 
they were irrelevant because y s < —2 for s > 1, one sees that the anomaly, 
beeing quadratic in s at first order in e and cubic in s at second order, could 
change this conclusion. Fixing d at physical dimension three, and say N = 3, 
we see that the s = 2 operator looks already relevant and it is a fortiori true 
for s > 2. Wegners's calculation is backed by the two-loop calculation (3) 
since the coefficient of e 2 is also positive. 

This is of course a very strange situation since it leads to the speculation 
that there would be an infinite number of relevant operators, whose coeffi- 
cients would have to be tuned before one reaches criticality. On the basis of 
similar calculations for generalized sigma models, Kravtsov et al. ||, were 
led to conclude that the one-parameter scaling theory of Anderson et al. [ITS 
could be invalid. 

For the Heisenberg model there is no experimental or numerical indication 
that one could have more than one relevant operator. Therefore one would 
conclude that, if these higher derivatives operators were truly relevant, al- 
though there is only one fixed point at zero temperature in dimension two, 
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one would have to find a new fixed point in dimension greater than two 
which would describe the critical properties. However the 1/N expansion, 
for instance, does not show any other fixed point, although it interpolates 
smoothly between the 4 — d and d — 2 expansions. So for large N at least, 
this scenario with some new fixed point, would run against all the existing 
evidence. 

Therefore let us return to the existing series they are obtained by 

a standard procedure in which e is a parameter which goes to zero first. The 
previous paradoxical conclusions are based however on fixing e and letting 
s grow. Clearly it involves an inversion of limits. To make this point more 
explicit, let us focus on the existing large s information. For large s we can 
write the expansion (|3|) as 

y s = d - 2s[l -~x- l -x 2 + 0(x 3 )] (11) 

in which we have defined the parameter 

es 

X = —2 < 12 > 



Let us consider the function 

f(x) = l- l -x- l -x 2 + 0(x 3 ) (13) 

For large s, one can speculate that for higher orders in e as well, the result 
for y s is of the form 

y. = l- 2sf(x) (14) 

It is quite possible that for the large s the result has this form and we are 
trying to prove it fl6 |. The sign of f(x), for fixed e, large s, i.e. for large 



x, is crucial. Near x = f(x) is positive, but its asymptotic expansion at 
the origin, extrapolated as it is, gives a negative result. However we know 
how dangerous it is to study the sign of a function for large x based on the 
knowledge of its asymptotic expansion at small x. Even if this expansion was 
convergent, it would presumably have singularities somewhere, and the small 
x expansion would give no information on the sign of the function beyond 
the closest singularity. Furthermore, although there are no rigorous results, 
it is more than likely that the expansion in powers of x is only asymptotic 
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and not convergent for any x. The most naive extrapolation of f(x) , given 
the existing data, would be to replace it by a [1,1] Pade approximant, i.e. by 



f( x ) = + 0(x 3 ) (15) 

We certainly do not claim that this is in any way a valid representation 
of f(x); however for large x, it is certainly as trustworthy as the small x- 
expansion Of course the Pade approximation (|i5|) is positive for large 

x, and the relevance of the large s-operators becomes even more doubtful. 

This discussion does not constitute a proof either in any way. We wanted 
simply to stress how unilikely is the picture with an infinite number of rel- 
evant operators, and how weak are the mathematical assumptions on which 
these speculations are based. 

This does not mean that the non-linear sigma models are always perfectly 
sound. Let us mention some real difficulties. Indeed consider an O(N)- 
invariant Landau potential V(<ft 2 ) which would develop a minimum away 
from zero below some temperature, in other words which leads to a first 
order transition. In the low temperature phase of the model, there are are 
still (N — 1) Goldstone bosons, since the analysis requires simply a broken 
continuous symmetry, but not necessarily a second order transition. In the 
low temperature phase of this model, we would write a low temperature 
expansion, which would show a fixed point at non-zero temperature above 
two dimensions. Clearly it would be wrong to interpret this fixed point as a 
critical temperature. 

The situation may even be a little more subtle. We know of cases in 
which the Landau potential gives a mean-field second order transition, but 
the fluctuations drive it to first order. This seems to happen if, for the same 
0(iV)-symmetry, the field belongs to the adjoint representation of the group, 
rather than to the vector one. It means that <fi is an iV x iV symmetric, 
traceless, matrix, which transforms under the rotation u as u T (j)uj. The Lan- 
dau potential involves two quartic invariants, namely (Tr0 2 ) 2 and Tr0 4 . The 
renormalization group at one-loop, in this space of two coupling constants, 
shows at fist order in 4 — d, a runaway solution to an unstable potential, 
indicative of a first order transition. However the analysis of the Goldstone 
mode of the 0(N)- symmetry broken down to 0(N — 1) depends uniquely 
upon the Lie algebra of these two groups, and not on the representation of 
the group to which belongs the order parameter. Thus, there again, we would 
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find a fixed point in the d — 2 expansion, whereas it is not at all clear that 
the model has a second order transition. 

Therefore we do not claim that the non-linear sigma models are to be 
trusted as other well established field theories such as QED. Our goal was 
simply to spell out that the speculations on their possible breakdown for 
the classical Heisenberg model, or for the one-parameter scaling theory of 
localization, are based on an unreasonable mathematical analysis. 
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